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ABSTRACT 


This paper describes a theoretical and experimental study 
of closely spaced open-ended waveguides radiating through a 
conducting ground plane. The waveguides are large enough in 
some cases to allow two orthogonally polarized propagating modes 
in each waveguide. 

A first-order analysis is presented which is based upon 
the method of moments and uses a single-mode approximation to 
the aperture field for each polarization. An improved first- 
order analysis is also presented which uses a higher order 
mode solution for the self admittances. Close agreement with 
experiment is demonstrated. 
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SUMMARY 

This report describes a theoretical and experimental study 
of closely spaced open-ended waveguides radiating through a 
conducting ground plane. The waveguides are large enough in 
some cases to allow two orthogonally polarized propagating modes 
in each waveguide . 

A first-order analysis is presented which is based upon the 
method of moments and uses a single-mode approximation to the 
aperture field for each polarization. An improved first-order 
analysis is also presented which uses a higher order mode solu- 
tion for the self admittances. Close agreement with experiment 
is demonstrated. 


INTRODUCTION 

The realization of large-scale microwave phased arrays has 
focussed new interest on the problems of mutual coupling between 
radiating apertures. A previous survey of developments in this 
field has been included in an earlier paper (refs. 1 and 2) . 

There are two distinct approaches to phased array problems. 
Work on infinite arrays of waveguides has progressed rapidly to 
a point at which several authors (refs. 3, 4, 5, and 6) have used 
higher order mode representations of the aperture field and have 
obtained solutions that are capable of explaining certain anoma- 
lous behavior (resulting in deep nulls in the radiation pattern 
as the array is scanned). Farrell and Kuhn (refs. 3 and 4) point 
out that this behavior cannot be predicted using a single-mode 
solution. The second approach to the array problem assumes the 
point of view that each terminal pair leading to an antenna ele- 
ment is a single port and is coupled to every other port in the 
array by means of the mutual admittances or the array scattering 
matrix. Since the controlling integral equations are linear, 
this approach can be made completely rigorous if the mutual cou- 
pling terms or scattering coefficients are evaluated in the 
presence of the whole array, and if the computation includes all 
of the higher order modes. A solution for the coupling between 
two col linear coupled slots has recently been published which 
does account for a large number of modes (ref. 1). There is 
strong justification, however, to seek simpler yet accurate 



mutual coupling formulas for aperture antennas, since this method 
can be used directly for small arrays or for edge effects in 
large arrays; can be readily compared with experiment, and can 
be used to study infinite arrays directly or after transformation 
using the Poisson summation formula (ref. 7). An example of the 
use of this type of analysis to study the performance of large 
scanned phased arrays is included in reference 8. 

The number of modes required to attain a specified accuracy 
or to describe a particular phenomenon is strongly dependent on 
the formulation of the problem. This statement is best supported 
by the well-known variational formulation results (refs. 9 and 10) 
for the rectangular waveguide. These one-mode results are far 
more accurate than the one-mode solution to the integrated form 
of Lewin ' s basic integro-dif ferential equation (ref. 2). 

Another interesting point raised by the work of Farrell 
and Kuhn, and a particularly relevant one to the development of 
realistic mutual coupling formulas, is the existence of cross- 
polarized modes at the array face. These fields are always 
present for arrays of waveguides with finite walls, but if the 
waveguides can support two propagating orthogonally polarized 
modes, then large cross-polarized signals are reflected back into 
each waveguide at certain directions of scan even though the 
array was driven with only one polarization. In a related devel- 
opment, Amitay et al . (ref. 11) have recently developed a match 
optimization procedure for finite phased arrays which requires 
a knowledge of the inter-element coupling coefficients for both 
polarizations. The authors compared their theoretical work with 
an experiment, and used measured values of the coupling coeffi- 
cients to perform the optimization. 

This paper presents a study of mutual coupling between two 
waveguides which radiate through a common ground plane. The wave- 
guide walls are parallel to one another, but the relative wave- 
guide position and orientation are arbitrary. Both components of 
polarization are considered in the coupling formulas, and the end 
result is a solution for the scattering matrix of the four-port 
network. The results are compared with experiments and with 
those of a more accurate theory in the special case of collinear 
slot coupling. 


FIRST-ORDER SOLUTIONS 
INCLUDING TWO ORTHOGONAL POLARIZATIONS 

In this simplified formulation each of two orthogonal polari- 
zations is represented by a single mode within each waveguide. 

The waveguide orientation is shown in Figure 1. The fields in 
free space are represented by the free-space dyadic Green's 
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Figure 1.- Coupled waveguide geometry 


function, and the fields in each waveguide are represented by 
the transverse field (ref. 12) expansion. Details are given in 
Appendix I. Equating the tangential magnetic fields in the 
apertures gives the resulting set of integral equations (ref. 7 
p. 100; ref. 13, p. 177). At waveguide No. 1: 

z x[ii(°)hi(xi,yi) + I 2 (o)h 2 (x 1 ,y 1 )] 


= z X 


/l( r l' r p)'P * E ( x p'Yp)] 


dSj 


(1 



At waveguide No. 2: 


z x [l 3 (o)h 3 (x 2 ,y 2 )+ I 4 (o)h 4 (x 2 ,y 2 )] 


where 


A 

= 2 


x E 

P=1 / 2 





( 2 ) 


i( r k ' r p) = y(i r k' r pD - y( r k P ) 

and 

r n = V ( x - x i) 2 + (y - y{) 2 

r 12 = V ( x 2 + D x • x l) 2 + (y 2 + D y • y l) 2 
r 21 - V( x 2 +D x- x l) 2 + (yj + ° y - Yi) 2 

r 22 = V ( X 2 - * 2 ) 2 + (^2 -y2) 2 


In these equations, the field subscripts 1 and 2 refer to the 
normal (H^q) and orthogonal (Hqi) polarizations in waveguide 
No. 1, and the subscripts 3 and 4 refer to the normal and orthog- 
onal polarizations in waveguide No. 2. The position subscripts 
Xp and yp are indicated in Figure 1. The integrals over S p are 
therefore integrals over apertures 1 and 2. 

There are a number of ways of proceeding to a solution com- 
mencing with Eqs. (1) and (2). Levis (ref. 10) multiplied these 
equations in scalar fashion by the aperture electric field and 
was able to show that the variational form which he obtained led 
to admittance parameters which are stationary with respect to 
the variations in the assumed fields. L^on et al . , (ref. 13) 
scalarly multiplied these equations by (z xh^) for each present 
on the left side of Eqs. (1) and (2) . The resulting equations 
are exact, but the solution was carried out by assuming that only 
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one mode was present in the waveguide aperture and that the 
waveguide-backed ground-plane slot had a total magnetic field 
equal to the incident magnetic field (as is the case for scatter- 
ing by a hole in an infinite screen) . It can be shown that if 
the solution were carried out rigorously, the value l 2 / v l com- 
puted using this second assumption with the second waveguide 
terminated in a matched load will be (l+y 33 /yg )/2 times the 
correct value. The analysis presented here begins with the 
procedure outlined by Lyon et al., but includes the additional 
orthogonally polarized component. Moreover, instead of separ- 
ating the equations using the magnetic field assumption the 
four equations are written directly in terms of their admittance 
matrix, and the scattering coefficients are obtained by perform- 
ing a matrix inversion and using the well-known relations between 
the admittance matrix and the scattering matrix of an N-port 
network (ref. 12, p. 13). In accordance with this outline, Eq.(l) 
is scalarly multiplied consecutively by z x h^ and z x h 2 , and 
Eq. (2) by z x h 3 and z x E 4 . The resulting equations are inte- 
grated over their respective apertures and reduced by orthogo- 
nality to the four equations below. This procedure is an example 
of the method of moments (ref. 14). 


V 



< 

X 2 

v J 

>= 2 / ds i Ms' < 

p= 1^ 2 *S 1 lJ 8 P 

1 P 

> 

Vx^i) 

^ J 



(4) 


(5) 




'' h 3 (x 2,y 2 ') 




> 


= p= 5 2 *( ds *^ ds p I | — (r 2 ,r p)*[ z x E (Xp,yp)] 


i 1 ^ 


2 p 


k h 4 ( X2,y 2 j 


(6) 


(7) 


After expanding the admittance dyad, writina the fields in terms 
of the modal functions, and performing the numerical operations, 
these equations may be written using the admittance matrix 
notation 


I = yv 


(8) 
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where the admittance matrix coefficient satisfies reciprocity 
Y mn = y nm and aynmetxy relations Y 1;l = Y 33 , y 22 = y 44 , y 32 = y 14 , 
y 34 = y±2‘ With these relations, only six parameters need be 
defined to determine all the others. 


y ll 

= y l 

c i' f a 

(o , a ; 

ro,b) 


y 12 

= y l 

C 1*B 

(o,a; 

ro,b) 


y 13 

= Y]_ 

C 1*A 

( D x 7 a;D y ' 

b) 

y 14 

= y x 

C 1*B 

(D x,a? D y , 

b) 

y 22 

= Y 2 

C 2 ¥ A 

(o,b; 

o, a) 


y 2 4 

= y 2 

C 2 ^A 

(D ,t 

y 

>;D x' 

a) 


The parameters ¥ B , y^, Y 2 r c^_, and C 2 are defined in 

Appendix II. 

The admittance parameters computed in this manner are not 
stationary with respect to the aperture field, although the self 
admittances are identical to those derived by Lewin (ref. 9) 
and, therefore, do represent the first-order stationary solution 
for an isolated waveguide. Though this self admittance is known 
to be accurate for most waveguide radiators, the slot conduc- 
tance it predicts for the square waveguides considered later in 
this paper is up to 10 percent less than is predicted by the 
more accurate twelve-mode solution (ref. 1). The error in power 
coupling which this introduces is discussed in the following 
section . 

Having thus defined all of the admittances, the scattering 
parameters may be obtained by matrix inversion using the equation 

S = (I - y) (I + y) _1 (10) 


where 


■ mn 


y mn^ y m 


These scattering parameters will be the items of interest for 
the remainder of this paper. 
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Before progressing to a discussion of the data, it is appro- 
priate to mention several conclusions evident from a considera- 
tion of the symmetry of 'P A and ¥ B . The function Tg = 0 whenever 
D y or D x or both are 0. Therefore, y \2 is identically zero, and 
y 14 (and, hence, S 12 and s 14^ is zero when the slots are arranged 
in either collinear or parallel fashion. In addition, it can 
also be shown that 'i' B (D x , a;Dy,b) = ¥ B (Dy ,b;D x , a) . This symmetry 
was used earlier to require that y-^ 4 = ¥ 22 ' anc ^ Yl2 = Y34 = *-* / 
but it also requires that if a = b, the coefficients S ^2 an< ^ S-^ 4 
have even symmetry in <j> about the angle <J> = 45 degrees. In 
addition, all parameters have even symmetry in (j> about <f> -= 0 , 

<j) = 90 degrees and, hence, about <f> = 180 degrees. For this 


reason, if a = b, only the parameters S^-, , Si2/ s 13' 


and 


S^ 4 need 


be computed in the restricted range 0 <. $ <_ 9 0 degrees to com- 
pletely specify all other coupling parameters at all orientations 


THEORETICAL AND EXPERIMENTAL RESULTS 

An experiment was conducted to compare with theoretical 
r suits. All data were taken with the slots mounted on a 6-ft 
by 6-ft ground plane with the equipment located behind the ground 
plane. Absorbing material was used to minimize the reflections 
from the walls of the enclosing room and, in some cases (Figures 
3A and 4) , the frequency was swept and the data averaged while 
making power-coupling measurements in order to minimize the 
error due to reflections. These measurements were repeatable 
to within about 0.5 dB. 

Since commercial dual-mode couplers were used to measure 
both polarizations of coupling at the received waveguide, a 
technique was based upon a previously published article (ref. 15) 
and required that an electric field probe be mounted perpendicu- 
lar to the ground plane in order to establish a known phase 
relation between the two slots when they are excited together. 
When <(> 0 , the probe is placed midway between the slots on the 

line joining their centers and the probe reading is null when 
the slots are excited and in phase. When <j> = 0 , the electric 
field component E on this center line is zero because of sym- 
metry, and so the probe is placed along a perpendicular bisector 
of the center to center line but away from this zero-field posi- 
tion. For this configuration, the probe reading is null when 
the slots are excited with equal out-of-phase signals. Other 
details of the phase-bridge calibration are included in the 
reference . 

The phase reference for the orthogonally polarized coupled 
signal was obtained by exciting that waveguide and setting the 
electrical path length equal to that of the waveguide receiving 
the transmitted polarization. The resulting ± 180-degree phase 
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ambiguity was resolved to ±360 degrees by correlation with the 
theoretical results. Coupled phases could have been measured 
by this technique to an accuracy of about ±4 degrees, if all 
external reflections were eliminated, the actual phase error in 
the presence of reflections was about ±10 degrees. 

Figure 2 shows the amplitude and phase of the measured 
coupling and compares this with results computed using both the 
first-order theory and a more accurate nine-mode solution 
(ref. 1) . The two computed values of power coupling are nearly 
the same for the lower frequencies, but as the frequency is 
increased, the first-order solution predicts somewhat less power 
coupling than either that observed or computed using the higher 
order mode solution. This deviation was greatest for the closest 
spaced slots for which the higher order mode coupling curve devi- 
ated quite visibly from a straight line. This deviation has been 
shown (ref. 2 ) to be due to the increased amplitude of higher 
order modes in the aperture field. The measured phase of the 
coupled signal shown in Figure 2B follows the nine-mode solution 
very closely. The first-order solution predicts the phase to 
an accuracy of about 10 degrees except for the slots which are 
very closely spaced. 

Figures 3 and 4 show the variation of power coupling and 
phase with frequency for square waveguides with several different 
values of <f>. These data are compared with experimental results 
over only a limited frequency range ( 8 . 5-9. 7 GHz) because of the 
band-pass range of the dual-mode couplers. The measured coupling 
was always within about 1.5 dB of the theoretical data and exhi- 
bited a nearly identical variation. In order to determine how 
much of this error was due to using the first-order self admit- 
tance, the theoretical results were recomputed using a higher 
order (12-mode) self-admittance. These results are shown by the 
dotted line in Figure 3. The coupling computed by this improved 
first-order theory agrees to within about 1 dB with the measured 
data, and this is nearly equal to the absolute accuracy of the 
measuring system. The phases of the coupled signals are shown 
in Figures 3B and 4. Unfortunately, the phase of S ^3 was so 
nearly constant between <f> = 45 and 90 degrees, as was the phase 
of S 14 for <j> between 30 and 6 0 degrees, that these phase changes 
with the angle <f> could not be detected experimentally. The theo- 
retical curves reflect the same characteristics. 

These comparisons show that this simple theoretical treat- 
ment can yield very good quantitative information about mutual 
coupling. Figure 5 shows the dependence of S-^, S^, and 
S ^4 with position angle <f> for various slot separations with 
a/A = b/A = 0.6. This curve was computed using the improved 
first-order solution. The coefficient S^i undergoes little change 
as these parameters are varied, and so is not plotted. Its value 
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is about 0.146 - j 0.020. As R c /A is changed in Figure 5A, the 
coupling amplitude | S 3^ 3 | at <J> = 0 degree varies much faster than 
than at <j> = 90 degrees. This is so because this amplitude varies 
asymptotically as 1/R^ when (j) = 0 and as 1/R C when <p = 90 degrees. 
Figure 5A shows, however, that for relatively close waveguide 
spacings the minimum coupling occurs at an angle between <j> = 0 
degree and <j> = 9 0 degrees, and as the spacing is increased this 
minimum moves closer to its asymptotic limit <f> = 0 degree. This 
coupling variation with <f> is therefore a near field effect and 
so any attempt to describe the scattering parameter S 33 as the 
product of an angle variable and a space variable would fail to 
describe it completely. 

The second important feature evident in Figure 5A is that 
the amplitude of coupling into the orthogonal port | S 3 ^ 4 [ is 
larger than [S 33 I for some angles of <p . Therefore, this term 
must be considered when the waveguide geometry is such that this 
mode can propagate. The coefficient IS 32 I is also shown in this 
figure and behaves essentially like | S 3 ^ 4 [ as (j) is varied, though 
it is about an order of magnitude smaller. Both | S 3^2 I and J 4 J 
are symmetric about (f> = 45 degrees. 

Figure 5B shows variation in phase of S32' S33, and S 14 as 
the waveguide angle <J> is varied. Of particular interest in these 
curves is the nearly constant phase change between curves of 
constant R Q /A , indicating the expected phase delay with distance. 
The phase change in S12 is about double that of S13 or S34 for 
the same change in R c /A , thus emphasizing that S^2 does not arise 
at the face of waveguide No. 1 itself but is due to the round 
trip coupling from waveguide No. 1 to waveguide No. 2 and then 
back again. 


CONCLUSION AND COMMENTS ABOUT APPROXIMATE SOLUTIONS 

This study has shown that the first-order or improved first- 
order mutual coupling approximations can be used to provide an 
accurate and convenient formulation for the four-port scattering 
parameters of two arbitrarily located waveguides radiating through 
a ground plane. This model has also shown that the cross-polarized 
component of the coupled signal is of the same order or magnitude 
as the coupled signal having the transmitted polarization for 
certain slot orientations. 

Should a higher order solution for the mutual coupling param- 
eters be required, some additional evidence uncovered in this 
study will help in simplifying the analysis by avoiding the need 
of solving the full-vector diffraction probelm. Instead of com- 
puting the four-port scattering matrix, S 13 was computed assuming 
that only the two ports (1) and (3) were coupled, and then S 14 
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Amplitude and phase of the waveguide scattering 
parameters using improved first-order theory. 
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is computed assuming that only ports (1) and (4) were coupled. 

The coupling between ports (2) and (3) or (4) was computed in 
the same manner using the following equation derived from Eq. (10) 
for these approximate cases: 



- 2 y 


mn 


(1 + 


) ( 1 + 

min 


nn 


) - Y 


2 

mn 


This approximation cannot be used to compute the second- 
order cross-polarized terms S 12 a ^d S 34 because Y\2 an< ^ Y 34 are 
zero. The other coupling parameters computed in this manner 
differ so little from those computed using the full four-port 
amittance matrix that they could not be distinguished in the 
accompanying curves. The same two-port model can be used to 
compute the reflection coefficients S rrm presented at each junc- 
tion, but here it must be remembered that each of the two-port 
calculations contains the self -admittance term and so one must 
subtract 


mm 


Isolated 


( 1 - y ) 
2 mm 

( 1 + y_ ) 

mm 


from the sum of the two Sj^ calculations computed using the 
two-port networks . 

These approximate results are valid because the second-order 
orthogonal mode coupling coefficients (S-j _ 2 and S 34 ) are very 
small. Since they are not a result of the coupling model used 
here, it seems likely that a similar approximation can be used 
in the derivation of coupling formulas if higher order mode 
solutions are desired. 

Concerning the need for higher order approximations , the 
data presented in this paper show that the formulation used here 
should suffice to compute the parameters of most phased arrays. 
Whether it will suffice to describe the anomalous behavior which 
has been recently reported must be determined by further study. 
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APPENDIX I 


Throughout this paper, two orthogonal modes (H 10 and Hqi) 
are used to describe the fields in each waveguide. The field 
expansion used is that of Marcuvity (ref. 9, Ch. 1-2). The 
subscript notation is 


Waveguide No. 1: 


Waveguide No. 2: 


Et. 


Ht. 


Et, 


Ht, 


E V- (z) e ■ (x, ,y, ) ; 
i=l , 2 1 

E I (z)h.(x ,y ) 
1 = 1 , 2 


E V. (z)e. (x„,y 9 ) ; 
1=3,4 1 1 ^ 


E 

i=3 , 4 


I.(z)h.(x 2 ,y 2 ) 


The H mode functions e^ are defined by: 


h . = z x e . and e . = z x v . 'F . 
l l l t l 


where 


”\ 

Y 1 


< 


7TX . 


sm 


> - - 2a J 
r ~ WF < 


’F, 


irx. 


sm 


and 


r -\ 

^2 


r ■N 

7ry n 


>= - 


sm 




2b 
7T a \ 


Try. 




sin 


J 
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The fields in free space are expanded in terms of the free- 
space Dyadic Greens function r° using the following relations: 


B(r) = j2ajye / r°(r,r' 

P = 1 / 2 o 


) • (z x E)dS 


where 


E (r) = 2 2 / VG(r,r ' ) x (z x E)dS ' 

p=l , 2 ''S P 


r°(r,r') = (U + -^ V7 ) G (r , r 1 ) 

k 

o 


and 


G ( r , r ' ) = 


-jk Q | r - r ' 


4 it r - r ' 


r - r 


= >4x - 


m2 


x') + (y- y ' ) + z 


• \ 2 , 2 


The admittance dyad is defined below: 


APPENDIX II 


The constants and C 2 are dimensionless and defined below. 
The characteristic admittances for the various model waveguides 
are also defined below. 


c-, = 


1 2 ir 2 (a/X) (b/X) 




C 2 “ „ _2 


2u^ (a/X) (b/X) ^ 1 _r^ 2 


'l = y 3 = V 1 - [A]' 

'2 = y 4 = Vi' [jfe]' 


The parameters and ¥ are also dimensionless: 




r a/2X 

r b/2\ 

-a/2X 

-b/2X 

= 

/ dx 

/ dy 

/ dx' 

/ dy 


J-a/2\ 

J -b/2\ 

■'-a/2X 

J-h/2\ 


. ~ 2 
Try' ttx o ~ 

cos — cos — G 
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Ill 


where 


and 



These basic integrals may be simplified somewhat by integrating 
them by parts and using trigonometric identities. After consi- 
derable manipulation, they may be written in the forms shown 
below: 



where 

2 

f A (T,a) = ^-[-u(t - a/2X) - v] (a - b/2A) G a 
2 

+ -j-[u(t - a/2A) + v] (a+ b/2A)Gg 
2 

+ — u ( x + a/2 A ) - v] (a - b/2A ) G c 

2 

+ -^“[u ( t + a/2 A ) + v] (a + b/2A) G^ 

and 

f B (T,a) = cos ^ cos E7X [ G A “ G B - G C + G d] 
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I 


Other parameters are defined below: 

G fl - G(D x/A ,T +a /2A;D y/v a + b/2X) 

g b - G (Vx' T+a/2X:D y/X'°‘ b/2X ) 

G C = G ( D x/> ' T ~ a / 2 K ' D y/X ' ° + 

G D - G(D x/> ,T-a/2A;D y/A ,a-b/2x) 

u = [4 - (A/a) 2 ] sin 
v = [4+ (A/a) 2 ] ^ oos ^ 
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